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Abstract 

Background. Acidification is one of the important impact cat¬ 
egories for life cycle impact assessment. Although its charac¬ 
terization has progressed during this decade through the em¬ 
ployment of midpoint approaches, only limited studies of 
endpoint approaches have been performed. 

Objective. This study aimed at developing damage function of 
acidification for terrestrial ecosystems in Japan. Damage func¬ 
tion expresses a quantitative relationship between the inventory 
and endpoint damage. 

Methods. The geographical boundary was limited in Japan both 
for emission and impact. In this study, sulfur dioxide (S0 2 ), ni¬ 
trogen monoxide (NO), nitrogen dioxide (N0 2 ) (NO and N0 2 
collectively mean NOJ, hydrogen chloride (HC1), and ammo¬ 
nia (NH 3 ) were considered as major causative substances of acidi¬ 
fication. Net primary production (NPP) of existing vegetation 
was adopted as an impact indicator of terrestrial ecosystems. 
The aluminum toxicity was adopted as the major factor of ef¬ 
fect on terrestrial ecosystems due to acidification. The leachate 
concentration of monomeric inorganic aluminum ions was se¬ 
lected to express the plant toxicity of aluminum. 

Results and Discussion. The results of damage function gave 
utilizable factors both for a midpoint approach and an endpoint 
approach; Atmospheric Deposition Factor (ADF) and Damage 
Factor (DF) applicable to the former and the latter, respectively. 
The ADF indicates an increase of H + deposition per unit area to 
an additional emission of causative substance. The additional 
emission corresponds to some alternatives in industry, not the 
baseline emission. The DF indicates the total NPP damage in all 
of Japan due to the additional emission of causative substances. 
The derived NPP damage is on the order of one millionth of the 
NPP itself. HC1 and NH 3 showed larger ADFs and DFs than 
that of S0 2 and NO x . The reason was ascribed to the relatively 
large source-receptor relationships (SRR) of HC1 and NFIj. 
Flowever, since the method applied to determine the SRR of 
HC1 and NH 3 has larger uncertainties than that of S0 2 and NO x , 
attention is needed to handle the difference. 

Conclusion. The damage function easily defines the concrete 
NPP damage due to an additional emission. The impact indica¬ 


tor, NPP, also has an advantage in its mass unit that is directly 
summable through the entire impact categories. Expansion of 
endpoints, such as in aquatic ecosystems, material degradation, 
human health, and biodiversity aspects of terrestrial ecosystems, 
is an important subject for future work. Further, uncertain analy¬ 
ses for major parameters will provide helpful information on 
the reliability of damage function. 


Abbreviations: ADF-Atmospheric Deposition Factor; AP-Acidi¬ 
fication Potential; Cl-Category Indicator; DF-Damage Factor; 
DI-Damage Indicator; EF-Effect Factor; LCA-Life Cycle Assess¬ 
ment; LCIA-Life Cycle Impact Assessment; LIME-Life-cycle 
Impact assessment Method based on Endpoint modeling; NDVI- 
Normalized Difference Vegetation Index; NNR-Non-Neutrali- 
zation Ratio; NPP-Net Primary Productivity; NSS-Non-Sea Salt; 
SRR-Source-Receptor Relationship 

Keywords: Acidification; acidifying substance; aluminum tox¬ 
icity; causative substance; damage function; Japan; life cycle im¬ 
pact assessment; net primary productivity; terrestrial ecosystems 


Introduction 

Acidification has been involved in many LCA studies as an 
impact category or an impact pathway, which is one of the 
important environmental issues. Acidification is caused by 
atmospheric deposition of acidifying substances such as 
sulfuric acid (H 2 S0 4 ) and nitric acid (HN0 3 ), mainly gener¬ 
ated from oxidized sulfur and nitrogen which is largely con¬ 
tributed from anthropogenic emitters. 

Early LCA studies based on midpoint approach employed a 
characterization factor, termed acidification potential (AP) 
as originally developed by Heijungs et al. (1992), to evalu¬ 
ate emissions of causative substances of acidification. AP 
only gives relative information on emission impact of causa¬ 
tive substances on the basis of SO z equivalent. Although 
several later studies (e.g. Potting et al. 1998, Huijbregts 1999) 
improved the circumstances considering regional sensitivity 
of ecosystems and atmospheric fate of causative substances, 
these were still confined to an equivalent load without in¬ 
formation on concrete damage caused by acidification. 
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Fig. 1 : LCIA framework of LIME 


Table 1: Entreatment of 'acidification': in several LCIA frameworks 


LCIA framework 

Causative substance 

Target endpoint relevant to terrestrial ecosystems 

Damage indicator 

Target region 

Eco-indicator 99 a 

SOx, NO x , NH 3 

Species extinction (vascular plants) 

PDF (species) 

Europe 

EPS Version 2000 b 

SO x , NO x , H 2 S, hf, hci, nh 3 

Crop growth 

Wood growth 

Soil acidification (base cation) 

Species extinction (life) 

Mass (kg) 

Mass (kg) 
Capacity (eq) 

NEX (-) 

World 

ExternE c 

Total acidity (from SO 2 , NO x ) 

Agriculture 

Forest and natural ecosystems 

WTP (Euro) 

WTP (Euro) 

Europe 

LIME d 

S0 2 , NO x , NH 3 , HCI 

Primary production of vegetation 

NPP (kg) 

Japan 


a Eco-indicator 99 (Goedekoop and Spriensma 2000) includes eutrophication aspects. 

b EPS Version 2000 (Steen 1999) originally deals with the endpoints in Table 1 as impact categories. Therefore, these include all relevant pathways, 
not limited to acidification. 

c ExternE (Holland et al. 1999) includes direct effects of air pollutants such as ozone, SO 2 and NO x . 
d LIME involves this study, and it deals with regional air pollution as another category. 


Given the circumstances, several LCIA frameworks tackled 
the acidification problem applying endpoint approaches, e.g. 
Eco-indicator 99 (Goedekoop and Spriensma 2000), EPS (Steen 
1999), and ExternE (Holland et al. 1999). In addition, LIME 
was developed in Japan as the result of the LCA National 
Project of Japan funded by METI/NEDO (Inaba and Itsubo 
2002). LIME is a comprehensive LCIA methodology with 
endpoint approaches, which also involved acidification in its 
framework (Fig. 1). Table 1 summarizes the features of each 
framework narrowed down to category endpoints of acidifi¬ 
cation regarding terrestrial ecosystems. It shows that the han¬ 
dling of causative substances, target endpoints, and damage 
indicators are largely different among the frameworks. 

This study, forming a part of LIME, focused on the impact on 
terrestrial ecosystems with priority among various endpoints 
regarding acidification. A reason was that terrestrial acidifica¬ 
tion tends to precede aquatic acidification, namely inland water 
is acidified after the depletion of acid neutralization capacity 
of its watershed. Moreover, since the major part of land is 
covered by vegetation, acidification has a potential to induce 


remarkable damage to terrestrial ecosystems, i.e. it is a pri¬ 
mary producer of ecosystems. As liming activities (e.g. 
Henrikson et al. 1995; Meiwes 1995) indicate, much work 
and time are necessary for recovering from acidified condi¬ 
tions. Hence, terrestrial ecosystems are likely to suffer the larg¬ 
est impact due to acidification among various endpoints. 

What is a suitable indicator for evaluating acidification im¬ 
pact on terrestrial ecosystems? LIME selected two safeguard 
subjects regarding ecosystems, i.e. primary production and 
biodiversity (see Fig. 1). Although both primary production 
and biodiversity are important safeguard subjects, biodi¬ 
versity aspects were excluded from this study. Biodiversity 
impact induced by acidification, however, was difficult to 
quantify due to the lack of information regarding inherent 
tolerance against acid for all relevant species and the risk of 
extinction due to acid stress. In this study, therefore, we only 
aimed at quantifying the effects on primary production. 
LIME adopted NPP, a rate of primary production utilizable 
in food webs, as an impact indicator for terrestrial ecosys¬ 
tems in relation to primary production. NPP is defined as 
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the difference between gross photosynthesized energy and 
respiratory-consumed energy by plants; it is expressed by 
mass as dry matter per unit area and time (e.g. kg ha -1 yr -1 ). 
It is considered that NPP is adequate as an impact indicator 
of terrestrial ecosystems, since NPP stands not only for the 
fundamental energy supply in food webs, but also for the 
indispensable oxygen supply in the biosphere. Meanwhile, 
the existing biomass of a primary producer, as a candidate 
for an indicator, denotes a stock energy only at a specific 
time point. Accordingly, existing biomass is inadequate to 
express the sustainability of terrestrial ecosystems. Although 
an LCIA study for land use (Lindeijer et al. 1998) employs 
NPP as an impact indicator, no similar study was found in 
the field of acidification. 

The purpose of this study is to develop the damage function 
of acidification for terrestrial ecosystems; damage function 
is defined as quantitative relationship between emission of 
causative substance and endpoint damage. The prototype 
of our study was shown in Hayashi et al. (2000). 

1 Methods 

1.1 Acidifying substance and its causative substance 

Major acidifying substances in the atmosphere originating 
from artificially emitted air pollutants are H 2 S0 4 and HN0 3 . 
However, the contribution of HC1, i.e. non-seasalt chloride 
(NSS-C1), to dry deposition was observed in Japan as well 
(Takahashi et al. 2001). Furthermore, NH 3 has a potential 
of univalent acid in charge balance corresponding to nitrifi¬ 
cation in soil, though NH 3 is a base in the atmosphere. 

Main causative substances of these acidifying substances are 
S0 2 , NO x (NO and NO z ), HC1, and NH 3 . A load of acid is 
expressed by a load of proton (H + ), as determined by the 
charge balance between total anion and total cation, exclud¬ 
ing H + on the basis of the Arrhenius acid-base theory ex¬ 
pressed with major ion species in the atmosphere by, 

= [ soa-] + [ no 3 -] + [ c/-] - [nh 4 + ] 

- [Na*] - [K-] - [Mg 1 *] - [Ca 1 *] (1) 

where each term has a unit of charge equivalent concentration 
(e.g. eq L _1 ), [H + ]° act is actual acidity without nitrification, and 
other terms correspond to each ion species, following the se¬ 
quence in the right side of Eq. 1, sulfate, nitrate, chloride, 
ammonium, sodium, potassium, magnesium, and calcium. 

As mentioned above, nitrification of deposited NH 4 + acts as 
an acid, with 2 mol of H+ per mol of NH 4 + , 

NHf + 2 °2 -> N Of + H 2 0 + 2H * (2) 

Inserting Eq. 2 into Eq. 1, actual acidity with complete ni¬ 
trification after deposition, [H + ] act , is expressed by, 

m act = [ SCV-] + [ NOf] + [ CM + [NHf] 

- [Na + ] - [K*] - [Mg 1 *] - [Ca 1 *] (3) 


In the context of acidification, atmospheric deposition loads 
of S0 4 2- , N0 3 ~, Ch, and NH 4 + contribute as sources of po¬ 
tential acidity. Therefore, these ionic species were consid¬ 
ered to be potential acidifying substances in this study. Fur¬ 
thermore, these acidifying substances were considered as only 
the source of acidity in this study. And other aspects like 
nitrogen as a nutrient were excluded from consideration. In 
addition, direct effects of air pollutants on terrestrial eco¬ 
systems such as oxidants and S0 2 itself were excluded from 
the impact category of acidification. As shown in Fig. 1, 
these effects were considered in the impact category of pho¬ 
tochemical oxidant and regional air pollution, respectively. 

1.2 Cause of damage 

There are many schools of thought on causes of forest de¬ 
cline (Hinrichsen 1986). In related to indirect effects caused 
by acidification, relevant causes are general-multiple stresses, 
aluminum toxicity, magnesium deficiency, excess nutrient 
such as nitrogen, heavy metal dissolution, and growth-alter¬ 
ing organic chemicals. Among these various causes, the plant 
toxicity of inorganic aluminum ion, in particular Al 3+ (Jenny 
1961), was selected as the primal factor of acidification im¬ 
pact on terrestrial ecosystems. It partly has to be ascribed to 
data availability in order to quantitatively relate causes and 
effects, though the aluminum toxicity is perhaps a serious 
factor. To express the aluminum toxicity, Al 3+ concentration 
in soil solution was adopted. 

1.3 Outline of damage function derivation and definition 
of factors 

Fig. 2 shows a schematic diagram of the damage function of 
acidification. The geographical boundary was limited to Ja¬ 
pan, namely both emission and impact only in Japan were 
considered. The damage function of acidification for terres¬ 
trial ecosystems consists of two factors. One is termed At¬ 
mospheric Deposition Factor (ADF) and reveals an average 



Fig. 2: Schematic diagram of damage function of acidification for terres¬ 
trial ecosystems 
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increase of H + load per unit area in Japan due to additional 
emission of causative substance. The other is termed Dam¬ 
age Factor (DF) which gives a total NPP decrease in all of 
Japan due to the additional emission of causative substances. 

ADF was determined using the fate of an emitted causative 
substance expressed as the product of Source-Receptor Rela¬ 
tionship (SRR) by Non-Neutralization Ratio (NRR) of acid 
in the atmosphere. SRR denotes the mean deposition fraction 
of an acidifying substance in Japan due to the emission of its 
causative substance in Japan. ADF is also utilizable for a mid¬ 
point approach; the summation of the products of ADF for 
each causative substance by its emission amount, namely, 
stands for the Category Indicator (Cl) of acidification. 

Meanwhile, DF was given by the product of ADF by Effect 
Factor (EF) that gives a total NPP decrease in all of Japan 
due to unit increase of atmospheric deposition load of H + 
per unit area, namely per ADF. EF was defined as the slope 
of the relationship between the increase of H + deposition 
load per unit area and total NPP decrease in the whole of 
Japan. The relationship was derived from quantitatively 
connecting the following processes; increase of H + load —» 
decrease of leachate pH —> increase of leachate Al 3 * con¬ 
centration —» decrease of plant growth rate —» NPP dam¬ 
age. DF is utilizable for an endpoint approach, namely the 
summation of the products of DF for each causative sub¬ 
stance by its emission inventory provides the Damage Indi¬ 
cator (DI) of acidification. DI denotes the concrete damage 
due to additional emissions of causative substances. 

1.4 Calculation of atmospheric deposition factor 
1.4.1 Source-receptor relationship of acidifying substance 

a. S0 2 and NO x 

In 1991, SRRs of S0 2 and NO x were derived from the esti¬ 
mation results of regional SRRs for 6 regions of Japan as 
the annual average values in Japan (Nakaminami 2000). 
These regional SRRs were calculated by modifying the 
Eulerian model for atmospheric transportation (Ikeda and 
Higashino 1997). The values of SRRs were estimated to be 
0.136 and 0.140 for SO z and NO x , respectively. 

b. HC1 

All of the NSS-Ch in the atmosphere was assumed to be HC1 
in this study. A value of SRR of HC1 in Japan was calculated 
using existing data on emission and deposition from different 
sources, since few studies on the SRR of HC1 were available. 

Waste incineration and fossil fuel combustion were consid¬ 
ered as HC1 sources, whereas the leakage from chemical in¬ 
dustry was neglected. The amount of generated HC1 at in¬ 
cineration was calculated by multiplying the amount of 
incinerated waste (SBJ 1998, 1999, 2000) by the average 
content of volatile chlorine in waste (NCUC-JWRF 1999) 
for each year from 1993 to 1997. Then, the amount of at¬ 
mospheric emission of incineration-origin HC1 was determined 
using the efficiency of flue-gas treatment facilities (Shigaki 
1998). As a result, the annual mean of HC1 emission from 
waste incineration was estimated to be 4.20 Gg-Cl ym 1 . Next, 
the amount of generated HC1 at combustion was calculated 
by multiplying the amount of oil and coal consumptions 


converted from the primary energy consumption of oil and 
coal (ANRE 1999) with the average Cl content in oil and 
coal (Okita 1982) for each year from 1993 to 1997. Then, 
the amount of atmospheric emission of combustion-origin 
HC1 was determined assuming the same removal efficiency 
of incineration-origin HC1. The annual means of HC1 emis¬ 
sion from oil and coal combustion were estimated to be 137 
and 3.16 Gg-Cl ym 1 , respectively. Hence, the annual mean of 
the total HC1 emission was estimated to be 145 Gg-Cl ym 1 . 

On the other hand, atmospheric deposition consists of wet 
and dry deposition. As the sum of wet and dry deposition, 
the annual deposition of NSS-C1 at the observatories of the 
Nationwide Survey of Acid Deposition was calculated for 
each year from 1993 to 1997 using the observation data of 
CF deposition (ADORC 2000) and the composition of ma¬ 
jor ion species in seawater (JWA 1990); the calculation was 
based on the assumption that all deposited Na + was of seasalt 
origin. NSS-C1 deposition was then calculated by subtract¬ 
ing the production of Na + deposition by the ratio of CF to 
Na* in seawater from CF deposition. These values, indicat¬ 
ing deposition at each observatory, were averaged for ur¬ 
banized and non-urbanized sites, and the total deposition in 
all of Japan was then calculated by summing the product of 
the averaged deposition by the area of each land classifica¬ 
tion. The annual mean of the total deposition of NSS-C1 
was estimated to be 57.0 Gg-Cl yr~F 

The SRR of HC1, derived by dividing the total deposition by 
the total emission, was 0.404. 

c. NH 3 

A value of SRR of NH 3 was calculated using existing data 
in this study similar to that of HC1. 

The estimated value of annual NH 3 emission in Japan in 
1994 was 460 Gg-N yr -1 (Kannari et al. 2001), including 
5% originating from natural sources. 

The annual mean of wet deposition of NH 4 + in Japan from 
1987 to 1989 was 336 Gg-N yr -1 (Fujita et al. 2000). This 
estimation was derived by using a grid system and corre¬ 
sponded to the area of 722,000 km 2 , excluding the area of the 
Nansei Islands, but including sea areas in each grid. There¬ 
fore, the value was simply separated for land using the land 
area in the same range, 370,000 km 2 , so that 172 Gg-N ym 1 
was then obtained. On the other hand, the annual mean of 
dry deposition of ammoniacal nitrogen at the observatories of 
the Nationwide Survey of Acid Deposition in 1994 was calcu¬ 
lated using the observation data (ADORC 2000). Then, the 
dry deposition in the whole of Japan was calculated in the 
same way for HC1. The annual mean of dry deposition of 
ammoniacal nitrogen was estimated to be 13.0 Gg-N yr 1 . 

The SRR of NH 3 derived from dividing the sum of wet and 
dry deposition by the total emission was 0.432. 

1.4.2 Non-neutralization ratio of acid 

Acidifying substances in the atmosphere are neutralized by 
coexisting basic substances. Here, it was assumed that acidi¬ 
fying substances which originated from additionally emit¬ 
ted substances were also neutralized as well as in the present 
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condition. Thus, NNR, defined as the ratio of actual acidity 
with nitrification to potential acidity in the atmosphere, was 
calculated as follows using the total deposition of H*, NH 4 % 
NSS-C1, NO,-, and NSS-S0 4 2 ~ at the observatories of the 
Nationwide Survey of Acid Deposition (ADORC 20000): 

NNR = ( H dep + 2 • NH 4dep ) / ( nssCl dep + NO idep 
+ nssS0 4dep + NH 4dep ) (4) 

where each term on the right side denotes the deposition of 
corresponding ionic species, following the sequence in Eq. 
4, proton, ammonium, non-seasalt chloride, nitrate, and non- 
seasalt sulfate. Each term in the right side has a unit of equiva¬ 
lent load (e.g. meq m -2 yr _1 ). Regarding the numerator, the 
reason why 2 NH 4 iep was added to H dep is that this H dcp was 
derived from the observed pH at which NH 4 * took a role as 
a base. Therefore, to determine actual acidity with nitrifica¬ 
tion, the generation of 2 H* per NH 4 + must be considered. 
In addition, NNR is principally related not to each ionic 
species directly, but to the acidity produced resulting from 
the total charge balance. 

NNR at each observatory was calculated for each year from 
1993 to 1997 and they were averaged for urbanized and non- 
urbanized sites. The weighted mean of NNR was determined 
using the respective area of urbanized and non-urbanized re¬ 
gions as the weight. Subsequently, the annual weighted mean 
of NNR from 1993 to 1997 was estimated to be 0.769. 

1.5 Calculation of effect factor 

1.5.1 Relationship between acidic deposition and soil 
acidification 

It was assumed that an increased H* load due to an addi¬ 
tional emission of acidifying substance completely contrib¬ 
utes to soil acidification, namely decreasing the value of soil 
pH. Based on the assumption, an increase of H* concentra¬ 
tion in leachate to an increased H + load was estimated; 
leachate is equal to vertically leaching water from soil. 
Leachate flux being used to determine an increase of H + 
concentration was simply expressed by, 

Q = Pr-ET (5) 

where Q, Pr, and ET denote annual leachate flux, annual 
precipitation, and annual evapotranspiration, respectively. 
Each term has a unit of mm yr- 1 . 


Q was calculated for each grid, termed 10 km grid, with the 
size of 5.0' for latitude and 7.5' for longitude. Pr was derived 
from the gridded dataset of meteorology as the normal values 
from 1961 to 1990 (MAJ 1996). The original dataset was 
also composed as a grid system with a ten-fold higher resolu¬ 
tion than that of the 10 km grid, and the arithmetic mean of 
Pr for each 10 km grid was calculated. On the other hand, ET 
was determined by applying the Thornthwaite method 
(Thornthwaite 1948), an empirical model as a function of 
monthly mean temperature and latitude. The monthly mean 
temperature in each 10 km grid was calculated, as well as the 
Pr employing the same dataset of Pr (MAJ 1996). 

An increase of H* concentration was determined by divid¬ 
ing an increase of H* deposition by Q, 

m, ne =(H, dd _ dep /Q)- 1(H (6) 

where [H*] mc and H add _ dep denote an increase in H + concentration 
(eq L _1 ) and increase in H* deposition (eq km -2 yr 1 ), respectively. 

Next, the representative soil type in each 10 km grid was 
determined as the most dominant soil type in the grid with 
the soil map of Japan at a scale of 1:1,000,000 (CSCN-GJP 
1990). The value of soil pH for each representative soil type 
was given as the median value of the corresponding soil type 
in the dataset of soil pH (Batjes 1995), data which was used 
to compile the soil unit of FAO-UNESCO Soil Map of the 
World (FAO-UNESCO 1990). 

Soil pH was assumed to be equal to leachate pH, and a value 
of leachate pH corresponding to an increase of H + concen¬ 
tration was then expressed by, 

P H,(lea) = - log ( 10 -P«*M + [H*] mc ) (7) 

where pH^lea) and pH 0 (lea) denote leachate pH with an in¬ 
crease of H* concentration and initial leachate pH, respectively. 

1.5.2 Relationship between soil acidification and Al 3+ 
concentration 

Although leachate pH affects leachate Al 3+ concentration, 
its quantitative relationship differs among soils. In this study, 
approximate models for each representative soil type were 
derived from existing data (So et al. 1999; Shioiri 1934; 
Kitagawa 1966; Yamaya 1968;Umemura 1968; NFRI1968; 
Maruyama 1995), which expresses the relationship between 
leachate pH and leachate Al 3 * concentration. The derived 
models are shown in Table 2. 


Table 2: Relationship between leachate pH and Al 3 * concentration derived for representative soil types in Japan 


Applied soil type in this study 

Exponential models for leachate pH - leachate Al 3 * concentration 0 

Soil or mineral type in original data 

Podzols, Podzolic soils b 

[A/ 3 *] 9.98E+3 x exp (-0.585 pH{/edj) 

Podzols 

Brown forest soils, 

Yellow-brown forest soils b 

[A/ 3 *] 9.13E+4 x exp (-2.19 pA(/ea)) 

Brown forest soils 

Red-yellow soils b 

[/4/ 3 *] 2.30E+6 x exp (-2.65 pH[/eO)) 

Red-yellow soils 

Andosols b 

[. A / 3 *] 5.33E+8 x exp (-3.69 pH[/eO )) 

Andosols 

Gley soils, Lowland soils, Lithosols b 

[A/ 3 *] 6.78E+2 x exp (-1.25 pH(lea)) 

Crystalline clay minerals c 


a phfjea ): leachate pH; [A/ 3 *]: leachate Al 3 * concentration (pmol L" 1 ). 

b These models were derived from existing data (So et al. 1999, Shioiri 1934, Kitagawa 1966, Yamaya 1968, Umemura 1968, NFRI 1968, 
Maruyama 1995). 

c The approximate equation of crystalline clay minerals was applied to gley soils, lowland soils, and lithosols, since no information on pH - Al 3 * 
relationships was obtained for these soil types. 
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1.5.3 Relationship between Al 3+ concentration and plant 
growth rate 

An experimental result with Pinus densiflora (Japanese red 
pine) (Lee et al. 1997) was applied to express the relation¬ 
ship between leachate Al 3+ concentration and growth rate 
of trees as dry matter. The derived model is as follows: 

GR = 8.00 • 1(H • [Al 3+ ] 2 - 0.0314 • [Al 3+ ] 

+ 100 ( [A/ 3+ ] <2000 ) (8) 

where GR and [Al 3+ ] denote relative growth rate of trees to 
control (%) and leachate Al 3+ concentration (pmol L" 1 ), re¬ 
spectively. 

1.5.4 NPP of existing vegetation and its damage 

The mean NPP (kg ha -1 yr -1 ) of each 10 km grid was calcu¬ 
lated using three map sources for existing vegetation, land 
use, and NDVI, and using a chart of NPP values for repre¬ 
sentative vegetation types. NDVI is an index of relative pho¬ 
tosynthesis activity. 

A grid system with 4.5" x 4.5" resolution, termed 100 m 
grid, was overlaid on a digital map of existing vegetation, 
the Dataset for GIS on Natural Environment (EAJ 1999), 
vector formatted, and the legend of the most dominant veg¬ 
etation in each grid was then extracted as the representative 
vegetation of the grid. The original legends of existing veg¬ 
etation of the map, approximately 900 legends, were inte¬ 
grated and simplified into 30 legends on the chart of NPP 
values (Iwaki 1981), and these NPP values were then allo¬ 
cated to each 100 m grid corresponding to its legend. The 
legend of each 100 m grid, however, was replaced by 'oth¬ 
ers', indicating highly artificial land, when the central point 
of the 100 m grid corresponded to the legend of a 'building 
site' or 'arterial road' of another digital map, the Detailed 
Gridded Map of Land Use in Japan with resolution of 3" 
for latitude and 4.5" for longitude (MLIT 2002). 

Furthermore, a digital map of NDVI, the Gridded NDVI 
Dataset with 36" x 36" resolution (CGER 2000), was over¬ 
laid on the produced NPP map of existing vegetation, and 
then the NDVI value at the central point of each 100 m grid 
was extracted. Subsequently, these extracted values were av¬ 
eraged for each of 30 legends. The ratio of the NDVI value 
at the central point of each 100 m grid to the averaged NDVI 
value corresponding to the legend at the point was applied 
to correct the fixed NPP value in the chart of Iwaki (1981). 
In addition, an NPP value of each 100 m grid was only cal¬ 
culated if the grid had all properties on existing vegetation, 
land use, and NDVI. 


Finally, the mean NPP (kg ha -1 ym 1 ) of each 10 km grid was 
determined by averaging NPP values for all the 100 m grids 
in the 10 km grid; here, the weight of 0.5 was applied to 
100 m grids which straddled two 10 km grids. 

NPP damage in each 10 km grid was calculated as follows: 

L-NPP act = NPP act ■ A ■ { GR 0 - GR 1 ) / GR 0 (9) 

where L-NPP act , NPP act , A, GR 0 , and GR 1 denote NPP dam¬ 
age in 10 km grid (kg yr- 1 ), average NPP of existing vegeta¬ 
tion in a 10 km grid (kg ha -1 ym 1 ), land area of 10 km grid 
(ha), growth rate at initial condition (%), and growth rate 
affected by additional H + load (%), respectively. 

The total NPP damage in the whole of Japan, L-NPP act (tot), 
was expressed by, 

n 

L-NPP ac , (tot) = X L-NPP act (,) (10) 

i 

where L-NPP act (i) denotes NPP damage of z’-th 10 km grid 
in Japan. 

2 Results 

2.1 Atmospheric deposition factor 

Due to an additional emission of its causative substance, a 
total deposition of acidifying substance in all of Japan was 
obtained by multiplying the additional emission of the sub¬ 
stance by its SRR. A mean deposition load per unit area was 
derived by dividing the total amount by the land area of Ja¬ 
pan. However, in the context of acidification, an additional 
emission amount of a causative substance was necessary to be 
converted into H + load. ADF (eq km -2 yjr 1 kg -1 ), which indi¬ 
cates the mean increase of H + deposition (eq km -2 yr- 1 ) in 
Japan by 1 kg yir 1 emission of causative substance, /, was 
expressed by, 

ADF (/ ) = 1 ■ 10 3 • SRR (j ) ■ MW ( j )~ 3 
■ VA {j ) • LA- 1 ■ NNR (11) 

where MW and VA denote molecular weight (g moH) and 
valence of acid of chemical species /, respectively. LA de¬ 
notes the land area of Japan (km 2 ). The concrete data to 
determine ADF are shown in Table 3, and the determined 
values of ADF are shown in Table 4. ADF is directly appli¬ 
cable to a midpoint approach. 


Table 3: Input data for determining Atmospheric Deposition Factors (ADFs) a 


Causative Substance 
(/) 

Source-receptor relationship 
(SRR) 

Molecular weight 
(MW) 

Valence of acid 
(VA) 

Land area of Japan 
(LA) 

Non-neutrali-zation ratio 
(NNR) 

S0 2 

0.136 

64.1 

2 

372,798 km 2 (1995) 

0.769 (1993-1997) 

Mean value 

NO x 

NO 

0.140 

30.0 

1 

N0 2 

46.0 

1 

HCI 

0.404 

36.5 

1 

nh 3 

0.432 

17.0 

1 


a ADF( j) [ eq kirf 2 yr 1 kg' 1 ] = 1 x 10 3 x SRR (j) x MW(JY 1 x VA (j) x LA' x NNR 
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2.2 Effect factor 

Connecting from Eq. 6 to Eq. 10 and inserting the following 
values of H add _ dep , 0, 1, 5, 10, 15, 20, 25, 50, 75, 100,150, 
and 200 eq km -2 yr -1 into Eq. 6, L-NPP act {tot) correspond¬ 
ing to each H add ^ ep was obtained. Plotting these data on the 
H adj-dep ~ L-NPP act ( tot ) coordinate system (Fig. 3), its slope, 
36859, was then adopted as EF, which indicates the total 
annual NPP damage in all of Japan corresponding to the 1 
eq km- 2 yr- 1 increase of mean H + load in Japan. 



Fig. 3: Relationship between increased proton deposition per unit area 
and total NPP damage in the whole of Japan 


2.3 Damage factor 

DF for each causative substance is shown in Table 4, which 
was obtained by multiplying ADF for each causative substance 
by EF. DF is directly applicable to an endpoint approach. 


2.4 Category indicator 

Cl indicates the total increase of H + load to the surface (eq krtr 
2 yr- 1 ) corresponding to an emission inventory. Cl is given by, 

S0 2 , NO, N0 2 , HCl, NH 3 

C/= £ Inv O’) ' ADF (j) (12) 

j 

where Inv { j ) denotes an emission of causative substance, j 
(kg yi -1 ). 

2.5 Damage indicator 

DI indicates the total NPP damage in all of Japan (kg-NPP 
yr -1 ) corresponding to an emission inventory. DI is given by, 

so 2 , no , no 2 , HCl, nh 2 

DI = X Inv(j)-DF(j) d3) 


3 Discussion 

3.1 Category demarcation 

There are some difficulties in clearly demarcating impact cat¬ 
egories regarding air pollutants but excluding chemicals, i.e. 
acidification, eutrophication, photochemical oxidant, and re¬ 
gional air pollution, since causative substances and endpoints 
of these categories are common and/or interactive. However, 
relevant processes of each impact category are largely different 
both in pathways and the seriousness of their impact; it forces a 
demarcation of impact categories regarding air pollution. 


Table 4: Damage function of acidification for terrestrial ecosystems 
Type 1: Emission unit is in kg of chemical species yr" 1 




Emission unit 

Midpoint approach 

Endpoint approach 

Causative substance 

ADF a 

AP b 

DF c 

so 2 

kg-S0 2 yr"' 

8.78 e-6 

1.00 

0.323 


NO 

kg-NO yr -1 

9.61 e-6 

1.09 

0.354 

NO* 

NOs 

kg-N0 2 yr"' 

6.27 e-6 

0.71 

0.231 

HCl 

kg-HCI yr' 1 

2.28 e-5 

2.60 

0.842 

nh 3 

kg-NH 3 yr"' 

5.24 e-5 

5.97 

1.93 


Type 2: Emission unit is in kg of S, N, or Cl yr’’ 




Emission unit 

Midpoint approach 

Endpoint approach 

Causative substance 

ADF® 

AP b 

DF c 

so 2 

kg-S yr-' 

1.75 e-5 

1.00 

0.646 


NO 

kg-N yr’ 

2.06 e-5 

1.18 

0.729 

INUx 

N0 2 

HCl 

kg-CI yr" 1 

2.35 e-5 

1.34 

0.866 

nh 3 

kg-N yr"’ 

6.36 e-5 

3.63 

2.34 


a ADF (Atmospheric Deposition Factor): Utilizable for Category Indicator (Cl) with emission inventory of substance / inv{J) [kg]. Cl= Z/nirJ x ADF 
(/). Cl has a unit of eq km" 2 yr" 1 which indicates a mean increase of IT deposition load due to additional emissions of causative substances. 
b AP (Acidification Potential): Equal to SO 2 equivalent ADF. 

c DF (Damage Factor): Utilizable for Damage Indicator (DI) with inv(j) [kg]. D/= XlnvlJ) x DF(/i. DF{J) = ADFfJ) x EF. DI has a unit of kg yr" 1 which 
indicates the total NPP loss in Japan due to additional emissions of causative substances._ 
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The consensus of handling these categories in LIME is as fol¬ 
lows: direct effects of air pollutants were divided into photo¬ 
chemical oxidant and regional air pollution considering strong 
acute toxicity of oxidants and chronic toxicity of other air pol¬ 
lutants, and then indirect effects of air pollutants were involved 
in acidification. Eutrophication through the atmospheric depo¬ 
sition was also a subject of acidification; eutriphication of closed 
water body was independently categorized (see Fig. 1). How¬ 
ever, eutrophication caused by atmospheric depositon was ex¬ 
cluded from quantification of damage at the present due to the 
difficulty in evaluation of damage, since eutrophication has as¬ 
pects of both merit and demerit to ecosystems. 

3.2 Cause of impact and impact indicator 

Aluminum toxicity is well known as one of the hypotheses 
of forest decline (Hinrichsen 1986) that was selected in this 
study as the main cause of indirect effects on terrestrial eco¬ 
systems induced by acidification. However, there is room 
for considering other causes of impact regarding indirect 
effects of soil acidification such as base cation deficiency 
and heavy metal toxicity. 

The concentration of leachate Al 3 * without threshold value, 
therefore being essentially toxic, was applied to express the 
aluminum toxicity in this study. Another index, e.g. the BC/A1 
ratio (Sverdrup and Warfvinge 1993), the molar ratio of base 
cations (K+Mg+Ca) to aluminum with threshold value, is 
ponderable for consideration in a future study. At the present, 
however, no spatial information on leachate concentration of 
base cation was available in all of Japan. It should be noted 
that only one tree species was involved to express the relation¬ 
ship of aluminum toxicity to plant growth, which is a weak 
point in the method. This is a subject to be improved. A dis¬ 
cussion of the representativeness of this relational expression 
of one tree species is to be seen in section 3.4. 

NPP was selected as the impact indicator of terrestrial eco¬ 
systems in this study. As shown in Fig. 1, there are two safe¬ 
guard subjects on ecosystems in the framework of LIME, 
i.e. primary production and biodiversity. NPP is an indica¬ 
tor expressing primary production. As mentioned in the in¬ 
troduction, biodiversity was excluded from the impact cat¬ 
egory of acidification, regretfully due mainly to a lack of 
information. Acidification may also possiblly affect 
biodiversity through its effects on interspecies competition, 
due to the differences in sensitivity against acidic condition. 
The perturbation of interspecies competition, however, is dif¬ 
ficult to model. Because insufficient information is available 
for all aspects, not only acid tolerance, but also competition 
mechanism among various species, and the chain or conse¬ 
quence of the change triggered by the first impact, it is very 
difficult to quantitatively evaluate the change in biodiversity 
quality, since it does not always entail extinction. 

We should discuss why the buffering capacity of soil to neu¬ 
tralize acid was excluded in this study. This property is never 
avoidable when dealing with inland water impact due to 
acidification, since inland water is what exudes through its 
watershed, i.e. soil and surface geology. In general, the total 
of these buffering capacities are large enough to neutralize 
substantial acid. This means that the low level of acid load 
which we consider in LCIA may bring about a negligible 


effect. Regarding terrestrial ecosystems, the impact of alu¬ 
minium toxicity on plant takes place on the root system hin¬ 
dering such features as root elongation and nutrient uptake 
(Cronan 1994). Fine roots, taking on the primary role of 
nutrient uptake, are mainly distributed in the surface layer 
of soil. At the surface layer, atmospherically loaded acid can 
affect plant roots before being neutralized. Because the sur¬ 
face layer is also the boundary between air and soil, and the 
input acid then has a chance to induce affects rapidly. Still, it 
is natural that some of the acid is neutralized even at the sur¬ 
face layer. However, it was difficult to mechanistically deter¬ 
mine how much of this acid will be neutralized. Thus, the 
buffering capacity of soil to acid was excluded in this study. In 
this sense, our estimation is related to the side of safety. 

3.3 Geographical boundary 

This study limited the geographic boundary of acidification 
within Japan to both the emission of causative substances 
and the impact on terrestrial ecosystems. In fact, acidifica¬ 
tion frequently becomes a transboundary issue typically 
emerged in Europe. Japan is an emitter of causative sub¬ 
stances of acidification and simultaneously a receptor of 
acidifying substances in East Asia. It is no doubt that atmos¬ 
pheric deposition in Japan is also effected by the pollutants 
from surrounding countries, e.g. China and South Korea. 
This could be involved in the baseline of acidic deposition if 
observation data were applied to determine this baseline. In 
this study, instead of directly including the baseline to esti¬ 
mate EF, however, the initial soil pH was assumed to be at 
equilibrium to the baseline of acidic deposition. Meanwhile, 
the aim of LIME was to evaluate additional effects caused 
by some alternatives in Japanese industries. Hence, the first 
purpose of this study was to evaluate effects due to addi¬ 
tional acid loads within Japan. In fact, however, there are 
many Japanese companies that have moved their produc¬ 
tion bases to abroad such as China. Thus, the transboundary, 
source-receptor relationships among countries in East Asia 
should be involved in the method in the future. 

3.4 Parameter representativeness 

SRR directly and proportionally affects the damage function 
of acidification so that it is important to know the fluctuation 
range of SRR for each causative substance, which changes 
year by year affected by conditions of meteorology and emis¬ 
sion. The SRRs determined for respective causative substances 
in this study were fixed values due to limited data availability; 
its representativeness and uncertainty are subjects to be ex¬ 
amined. The SRRs of HC1 and NH 3 were three-fold larger 
than that of S0 2 and NO x . The estimation method applied to 
determine the SRRs of HC1 and NH 3 was relatively rough. 
Thus, the uncertainties of their values are larger than that of 
S0 2 and NO x , which brings about difficulties in comparison. 
However, the difference can be attributed, in part, to the prop¬ 
erties of HC1 and NH 3 , i.e. highly water soluble and then high 
velocities of dry deposition (Puxbaum and Gregori 1998). 
Although S0 4 2 - and N0 3 “ generated from oxidizing S0 2 and 
NO x are also highly water soluble, oxidization needs time to 
progress. Therefore, S0 2 and NO x tend to be transported long- 
range during this time, which results in smaller SRRs. 
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Although Japanese red pine is one of the common tree species 
in Japan, it is rather a pioneer species and its forests cover 
small areas in Japan. However, no similar model to Eq. 8 was 
derived regarding other tree species that consist of dominant 
vegetation in Japan. Accordingly, regardless of vegetation leg¬ 
ends, Eq. 8 was applied to express the relationship between 
leachate Al 3+ concentration and the growth rate of a tree. The 
sensitivity of Japanese red pine to aluminum tended to be higher 
than that of Fugus crenata (Japanese beech) from experiments 
on the BC/A1 ratio (Izuta et al. 2001). Thus, it is possible to 
overestimate NPP damage in Eq. 8. 

The correction using NDVI for NPP values (see 1.5.4) was 
to describe the difference of NPP by climate and landform, 
since the original NPP dataset (Iwaki 1981) only gives fixed 
values. This correction method was developed and discussed 
in Nakagawa et al. (2002) as a part of LIME. The result 
could reflect the changes in NPP by climate, though it per¬ 
haps has some uncertainty. 

3.5 Results of damage function 

The DFs obtained in this study showed the similar order of 
values compared to emission of causative substances on mass 
basis (see Table 4). It was just a coincidence coming down 
to the size of geographical boundary, limited to Japan. 
Namely, a value of EF becomes large if the range of geo¬ 
graphical boundary was extended; it directly results in an 
increase of land area where its NPP damage should be 
summed. On the other hand, a value of ADF is relatively 
stable regardless of the size of geographical boundary, un¬ 
less a value of SRR changes remarkably. 

The relationship between increased H + deposition and total 
NPP damage in Japan, i.e. EF, showed complete linearity (see 
Fig. 3). However, even the uppermost value of increased H + 
deposition, 200 eq km -2 yr -1 , accounts for only 0.2% of the 
baseline H + load in Japan, around 100,000 eq km -2 ym 1 (1000 
eq ha -1 yr -1 ). Namely, the range of increased H + deposition is 
within very low loads and the relationship could keep its lin¬ 
earity, although non-linear relationships were used for deter¬ 
mining the EF such as H* concentration and pH, pH and 
aluminum concentration, and aluminum concentration and 
plant growth. It is, however, important to notice that the low 
H + load, 200 eq km' 2 yr -1 , corresponds to 74 Meq yr 4 of 
additional H + load in the whole of Japan. It is equivalent to 

3.5 Gg of S0 4 2 ' or 4.6 Gg of N0 3 ' in charge balance, all to¬ 
gether no small amount. The usage of the method is mainly to 
estimate additional impact due to some alternation in indus¬ 
try. Accordingly, the range of increased H + deposition seems 
to be sufficient. More attention is also needed to evaluate the 
magnitude of total NPP damage. As shown in Fig. 3, the or¬ 
der of total NPP damage is 10 s kg yr -1 , which indicates the 
nationwide sum of NPP damage in Japan. This damage is 
equivalent to 3 kg km -2 yr -1 (0.03 kg ha' 1 ym 1 ), considering 
the area of Japan. On the other hand, NPP itself is on the 
order of 10,000 kg ha -1 yr" 1 . Consequently, the NPP damage 
which is involved in this study contributes to one millionth of 
the NPP itself, which sounds like quite a risk. 

NH 3 showed the largest DF and HC1 was next, although SO z 
and NO x are major causative substances of acidification. The 
result was ascribed to the relatively large SRR and low mo¬ 


lecular weight for NH 3 and the relatively large SRR for HC1 
(see Table 3). The result is likely to include the fact, but the 
uncertainties of DFs for NH 3 and HC1 are essentially larger 
than that of SO z and NO x . A lot of attention is needed to 
handle the difference of DFs, and the reliability of relevant 
parameters, SRR in particular, should be enhanced. In addi¬ 
tion, S0 2 and NO x perhaps dominate emission inventory in 
industry, since these substances are evidently generated in quan¬ 
tity by energy consumption compared to HC1 and NH 3 . 

3.6 Comparison with another study 

As shown in Table 1, it was difficult to compare our dam¬ 
age function directly to other endpoint approaches because 
of the differences in frameworks. Here, a comparison on 
the midpoint level between our result and the result of Hei- 
jungs et al. (1992) updated with Hauschild and Wenzel 
(1998) was carried out. The value of AP in this study for 
each causative substance, AP(j), was derived as follows: 

AP(j) = ADF(j)/ADF(S0 2 ) (14) 

Table 5 shows the comparison result. APs of S0 2 and NO x 
were similar in both study results. APs of HC1 and NH 3 , on 
the other hand, were different and our result showed larger 
values. The reason was ascribed to the SRR. The similar value 
of SRR for NO x to that of S0 2 (see Table 3) brought about the 
similar APs between both study results, although Heijungs et 
al. (1992) took no consideration of SRR. However, the larger 
values of SRR of HC1 and NH 3 than that of S0 2 (see Table 3) 
resulted in the larger values of AP (Table 5) in this study. The 
SRR is indispensable to determine AP, though it perhaps has 
considerable uncertainty. Simply because acidification impact 
takes place on the surface and not in the atmosphere. Accord¬ 
ingly, an improvement of the reliability of SRR is important. 


Table 5: Comparison of Acidification Potential (AP) 


Causative substance 

i AP" 

This study 

Heijungs et al. b (1992) 

so 2 

1.00 

1.00 

NO, 

NO 

1.09 

1.07 

no 2 

0.71 

0.70 

HCI 

2.60 

0.88 

NH 3 

5.97 

1.88 


a AP: kg-S0 2 equivalent kg' 1 . 

b The result of Heijungs et al. (1992): updated with Hauschild and 
Wenzel (1998)._ 


4 Conclusion and Future Work 

The damage function of acidification for terrestrial ecosys¬ 
tems easily gives concrete NPP damage due to an additional 
emission of its causative substance. The impact indicator, 
NPP, also has an advantage in its unit, mass, that is directly 
summable through the entire impact categories, and pro¬ 
vides a chance of comparison between other impact indices 
such as crop production and fishery production. 

Expansion of endpoints, not only NPP damage of terrestrial 
ecosystems but also reproduction damage of aquatic eco¬ 
systems, material degradation, human health with regard to 
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drinking water toxicity, and of course the aspects of 
biodiversity, is an important subject for future work to com¬ 
prehend the present LCIA framework of the impact category 
of acidification. Further, uncertain analyses for major rel¬ 
evant parameters will provide helpful information on the 
reliability of damage function. The SRR, in particular, plays 
an important role in determining the effects of a load of acid 
on the surface through the atmosphere. 
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